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In the parity-time-symmetric (PT-symmetric) Hamiltonian theory, the optimal 
evolution time can be reduced drastically and can even be zero. In this letter, 
we report our experimental simulation of the fast evolution of a PT-symmetric 
Hamiltonian in a nuclear magnetic resonance quantum system. The 
experimental results demonstrate that the PT-symmetric Hamiltonian system 
can indeed evolve much faster than the quantum system, and the evolution 
time can be arbitrarily close to zero. 
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1. Introduction 

The brachistochrone problem, i.e., the minimum time evolution between 
two states, is an interesting and important problem. In quantum mechanics, 
the brachistochrone between two states is bounded by the maximum difference 
of the eigenvalues of the Hamiltonian |TJ[21[31I11[S1[S1[71[H]- Brachistochrone 
has important applications. For instance, the time-optimal approach to the 
quantum algorithmic complexity has attracted much interest recently I10|. 

The Hermiticity requirement of a Hamiltonian guarantees that its eigenval- 
ues are real. It also implies that the evolutionary operator e~'^^* is unitary. 
However, Hermiticity is a sufficient condition but not necessary for real eigen- 
values, and the entire spectrum of a wide class of non-Hermitian Hamiltonians 
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can also be real. Among these Hamiltonians [TT] is a class that is PT symmet- 
ric. The PT-symmetric Hamiltonian has been investigated intensively in recent 
years, both in theory [HlIIllIlIiailSiniliailllinilSnilllllSlllllMlllSl 
[Ml [571 [Ml 123 [33 E] and in experiments [Ml [31 [31 [Ml [3S| ■ PT symmetries 
have been experimentally observed in table-top optical systems [321 [331 [32 155] 
and in spin-polarized Rb atoms [36]. 

The novel character of PT-symmetric Hamiltonians brings about many new 
features and may lead to interesting applications. Faster than Hermitian quan- 
tum mechanics evolution is one of these important aspects [12]. In this article, 
we design and carry out an experiment that simulates the time evolution of a 
PT-symmetric Hamiltonian with a nuclear magnetic resonance (NMR) quan- 
tum system. We build a system in Hilbert space that admits both unitary 
and non-unitary evolution, and observe the time evolution of a PT-symmetric 
Hamiltonian The experimental result shows that the minimal evolutionary time 
in a PT-symmetric system is faster than that in the Hermitian case, and can be 
arbitrarily close to zero. 
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2. Theoretical Frame 

A simple PT-symmetric non-Hermitian Hamiltonian for a two-level system 

According to Ref. [T^ . the largest and smallest eigenvalues are E± = 2s cos a 
and 0, respectively. The difference between them is 

Lu = E+ - =2s cos a. (2) 

Under the influence of e~R^*, the PT-symmetric system that is initially in 
|0) = (1 0)"^ will go to [12J 
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to evolve to state |1) = (0 1)"^. When a — 7r/2, it approaches zero, which is 
an impossible task for a regular Hermitian Hamiltonian. 

For comparison, the equivalent Hermitian Hamiltonian, Hq, was calculated 

m 

rr f s COS a s COS a \ ,^s, 
^0 = , (5) 

\^ s cos a s cos a J ^ 

the eigenvalues of which are E± = 2s cos a and 0, respectively. Here, E^ — E^ = 
w, which is exactly the same as that in the PT-symmetric case. 
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Figure 1: Quantum circuit for a PT-symmetric system. Operations from left 
to right are a single-qubit operation V ^ Co-c/i (a O-controUed C/i), C\-U2 (a 
l-controUed U-i), and a Hadamard operation H. When the ancillary qubit is in 
|0)a, the final state of work qubit is e^'s^^*|0)e- 



The evolution under this equivalent Hermitian Hamiltonian is given by 

e-i^«* M = e^^*"™*^" cos(^^t) \ , . 



J \ sin(^^t 

and the time it takes to evolve to the final state |1) is 



TO = (7) 



which is constant for a fixed uj. 



3. Construction of a PT-symmetric Hamiltonian System 

We now construct a system with a PT-symmetric Hamiltonian Eq. ([1]) and 
simulate its time evolution. The vital part of the simulation is achieved using the 
idea of duality quantum computing |371 138) . Duality quantum computing can 
be achieved by using an ancilla qubit using a conventional quantum computer 
|371 138| . The principle to simulate a duality quantum computing is shown in 
a quantum circuit in Fig. [1] Simulating the non-unitary evolution of a PT- 
symmetric quantum system requires the determination of the explicit forms of 
the operators in Fig. [TJ 

For a conventional quantum computer, the idea to use an extended space 
consisted of a system and an ancilla for simulating a non-unitary evolution in a 
PT-symmetric system by a unitary evolution in an extended space was proposed 
in Ref. f39]. This is similar to the our scheme used in the experiment in this 
work. 

The system we use contains two qubits: the work qubit e and the ancillary 
qubit a. A qubit is a two-level quantum system that is the building block in 
quantum information processing. The input 2-qubit state on the leftmost is 
|0)a|0)e. We then perform the V unitary operation, 

Y — f '^^^ 't'v — sin (j)v 
y sm(j)v cos(f>v 
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on the ancillary qubit, where 0y is 



arccos 
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Next, we apply two controlled unitary operations, 



where 
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and U2 ~ (^z- Finally, a Hadamard operation is utilized on the work qubit e. 
Here, t is the evolution time in the PT-symmetric system, s is a parameter in 
the Hamiltonian in Eq. ([T]) , and tu is the difference between the two eigenvalues 
of the Hamiltonian. 

After performing operations shown in Fig. [1] the final state becomes 



q 

V2 



1 



|0)ae-^^*|0)e+ |l)a-(cOS0vC/i-sin,^yC/2)|O), 



where q is 



^irts COS a 

e " 




, ujt 

'2h 



sin ^ ^ ^ 
2s sin a ^ 



(12) 



(13) 



which is a non-zero number and tends to l/\/3 as a — > — 7r/2 and t = t. If we 
observe the work qubit conditioned on the ancillary qubit to be in state |0)a, 
the evolution associated with the work qubit becomes 



|0)e, 



(14) 



which is the PT-symmetric Hamiltonian evolution. 

It is worth explaining the symbols. We use t to denote the evolution time in 
the PT-symmetric quantum system. The time it takes to complete the evolution 
in the work-ancilla two qubits system is designated as t. The time it takes for 
the PT-symmetric quantum system to evolve from |0) to |1) is indicated as r 
and the corresponding time it takes to complete the operations in the work- 
ancilla two qubits system represented by f. As we construct the PT-symmetric 
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system with a one-qubit-subspace of a two-qubit-Hilbert space, in which the 
sub-system evolves non-unitarily while the whole system evolves unitarily, the 
evolving time i of the whole system depends on the evolution time t of the PT- 
symmetric system, and vice versa. The evolving time r for the PT-symmetric 
system can approach zero when a tending to — -j, which is faster than the 
counterpart Hermitian system. 

4. Experimental Realization 

We simulated the evolution process in an NMR quantum system, C^'^-labeled 
chloroform that consists of two qubits. The C^'^ nuclear spin works as the work 
qubit and the proton nuclear spin works as the ancillary qubit. We begin from 
the state |0)a|0)e. We then evolve the PT-symmetric system to some time t by 
applying the corresponding operations given in Fig. [TJ Next, we measure the 
state of the work qubit conditioned on the ancillary qubit a being in state |0)a. 
By varying the instant we observe the state evolution of the PT-symmetric 
system. 

The following NMR notations are adopted. The free evolution of the two- 
qubit system for a period of X is denoted as [X] : 

[X]=e-'^''>^', (15) 

where J 215.23 Hz is the coupling constant between C^'^ and H^. A rotation 
of spin m through angle cj) about axis j is denoted as [4>\j^ ^ and [(/)]™ — e~*'^'^J 

The spatial-averaging method was used to prepare the pseudo-pure state 
|0)a|0)e. The single-qubit rotation operation V is realized by a pulse 

[20y]^ (16) 

on the ancillary qubit. Here, Co-(7i and C1-U2 ^''^ realized by the following 
two pulse sequences: 

^ ^ [f ]1, ^ (17) 



and 



TT TT TT TT TT 

-> [3]? ^ [2]% ^ [3]^ ^ [3]" ^ ^2^ 
respectively. Finally, the pulse sequence 



(18) 



T^Ty -> N"-., (19) 
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Figure 2: Typical spectra of the work qubit with a = —1:1%: (a) pseudo-pure 
state at the beginning of evolution; (b) final state after evolving for a time of r. 



is applied to implement the Hadamard operation on the ancillary qubit. The 
evolution is observed by looking at both the state of the ancillary and the work 
qubit. If the ancillary qubit is in state |0)a, the state of the work qubit gives 
the evolution under the PT-symmetric Hamiltonian. 

Because we are interested in the behavior of the system near a = — 7r/2, we 
restricted the parameter in the range of a e (— 0]. We performed a series of 
experiments with values of a at 0, — — ~^lf"- IH^) 

shows the spectrum of the work qubit in the pseudo-pure state. A right single 
upward peak in the spectrum indicates that the work qubit and the ancillary 
qubit are all in state |0). Fig. Elb) shows the spectrum of the work qubit for 
a = — ^ . The downward peak on the right indicates that the two-qubit state is 
|0)a|l)e, whereas the peak on the left corresponds to the ancillary qubit in state 
which is not of interest to the present study. 

We now examine the simulation when f = 0. Here, V operation becomes 
an identity operator and it does not require time to simulate. The two free 
evolutions [(/)(7i/(27r J)] pulses and the last single qubit pulse [— hi the 
Co-(7i in Eq. (|17p also do not require any time. The other pulses in Cq-jj^ and 
the whole pulse sequence of C\-u^ and H still require a constant (with respect 
to €) time to complete. Thus, the simulation performed in the two-qubit system 
still requires time t to complete, even though t is zero. 

Quantitative results for the evolution in the PT-symmetric system were ob- 
tained by performing quantum state tomography. Fig. [3] gives the density 
matrices of work qubit for a — — 0.48447r at the beginning it = 0), middle 
{t = t/2), and end {t — r) of the evolution. Fig. [3]Ja) is the density matrix at 
the beginning where the state is |0); Fig. [3l[b) shows the state in the middle 
of the process; and Fig. Wic) gives the final state. For comparison, we drew 
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Figure 3: State tomography of the work qubit with a = — 3l7r/64 = — 0.48447r: 
(a) pseudo-pure state |0)e, (b) middle state at t — t/2, and (c) final state. In 
each picture, the left side provides the experimental results, whereas the right 
side provides the theoretical results. (Online version in colour.) 
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Figure 4: Total experimental time f for simulating evolution from |0) to 
The red circles denote the PT-symmetric case, while the green squares indicate 
the Hermitian case. The connected lines are used as visual guides only. (Online 
version in colour.) 

the corresponding theoretical density matrices on the right side of each picture. 
Clearly, the experiments agree with theory very well. 

The total experimental time f to finish the simulation of the evolution from 
|0)e to |l)e in the two-qubit system is shown in Fig. |31 As a approaches — 7r/2, 
f decreases remarkably; however, it does not reach zero. As the PT-symmetric 
Hamiltonian is realized in a larger quantum system, f depends, not only on the 
evolution time t in the PT-symmetric system, but also on the time it takes to 
set up the PT-symmetric system. The relation between t in the PT-symmetric 
system and t in the two-qubit system is determined through the four operations 
shown in Fig. [T] Of the four operations, V and Cq-ui are dependent on t, 
whereas Ci-U2 and H are constant operations that are independent of t. In 
the simulation, t represents a parameter in determining the operations of V and 
Co~ui- 

The relation between r and a predicted in Ref. [12J appears in the data 
remarkably well. An evolution faster than the Hermitian Hamiltonian system 
evolution is clearly observed. The evolution time r taken by the PT-symmetric 
system to go from |0) to state |1) is shown in Fig. [5] In Fig. O the evolu- 
tionary time for the PT-symmetric Hamiltonian system approaches zero when 
a approaches — 7r/2. 

For comparison, we simulated the evolution of the equivalent Hermitian 
Hamiltonian Hq. The quantum circuit is similar to that for the PT-symmetric 
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Figure 5: Evolution time t versus a. The red circles indicate the PT-symmetric 
case, while the green squares denote the Hermitian case. (Online version in 
colour.) 



case shown in Fig. [T] However, we substituted Ui, U2 and H by 

r~r TT ( cos (^-£2^0 -isin(^^t) \ 

= = ( sin {^t) cos (i^^) ^ ) (20) 

and V"^ . respectively, where V'^ is the Hermitian conjugate of the V operator 
in Eq. ([H]). The total time f to implement the simulation and the evolution 
time r for the equivalent Hermitian system are obtained and shown in Figs. [4] 
and m respectively. The evolutionary time r for the Hermitian case is clearly 
constant regardless of the value of a. The faster-than-Hermitian evolution of 
the PT-symmetric system is evident. 

5. Conclusions 

In conclusion, we experimentally simulated the evolution of a PT-symmetric 
system in an NMR quantum system with two qubits. The faster than 
Hermitian quantum mechanics evolution of a PT-symmetric system predicted 
in Ref. [12] is clearly observed. When the parameter a approaches — 7r/2, the 
evolution time also approaches zero. 

When the difference between the large and small eigenvalues of a Hermitian 
two-level quantum system is fixed, the fastest evolving time is invariant for a spin 
flipping in Hermitian quantum mechanics. However, for PT-symmetric quantum 
system, the brachistochrone time can be varied by changing the parameters in 
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the Hamiltonian. It can not only accelerate the evolution as predicted in Ref. 
[12] and demonstrated in this work, but also decelerate the evolving time as 
shown in Ref. |41j . 

This work was supported by the National Natural Science Foundation of China 
(Grant No. 11175094) and the National Basic Research Program of China 
(2009CB929402, 2011CB9216002). 
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